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Small Scale Structure at High Redshift: III. The dumpiness of 
the Intergalactic Medium on Sub-kpc Scales ^ 

Michael Rauch^, Wallace L.W. Sargent'^, Thomas A. Barlow^, Robert F. Carswell'^ 

ABSTRACT 



Spectra obtained with the Keck HIRES instrument of the Lyman a forests 
in the lines of sight to the A and C components of the gravitationally lensed 
QSO, Q1422+231, were used to investigate the structure of the intergalactic 
t^^ ■ medium at mean redshift < z >~ 3.3 on sub-kpc scales. We measured the 

cross-correlation amplitude between the two Lya forests for a mean transverse 
separation of 120/i^Q^pc, and computed the RMS column density and velocity 
differences between individual absorption systems seen in both lines of sight. 
The RMS differences between the velocity centroids of the Lyman a forest lines 
[■^^ ■ were found to be less than about 400 ms , for unsaturated HI absorption lines 

Q . with column densities in the range 12 < logN(HI) < 14.13. The rate of en- 

ergy transfer into the low density IGM on a typical scale of 100 pc seems to be 
^ . lower by 3-4 orders of magnitude than the rate measured earlier for strong CIV 

I . metal absorption systems. The tight correlation between HI column density and 

^ . baryonic density in the intergalactic medium was used to obtain a conservative 

c/3 . upper limit on the RMS fluctuations of the baryonic density field at < ^ >= 3.26, 

namely, a /((A log p) ) < 3.1 x 10^^ on a scale of 110 h'^^'pc. The fraction of the 
Lj absorption lines that are different across the lines of sight was used to determine 

Vh ' the filling factor of the universe for gas which has suffered recent hydrodynamic 

disturbances. We thereby derived upper limits on the filling factor of galactic 
outflows at high redshift. Short-lived, short-range ancient winds are essentially 
unconstrained by this method but strong winds blowing for a substantial fraction 
of a Hubble time (at z = 3.3) appear to fill less than 20% of the volume of the 
universe. 
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generous support of the W.M. Keck Foundation. 
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^Institute of Astronomy, Cambridge University, Madingley Road, Cambridge CB30IIA, UK 
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Subject headings: intergalactic medium — cosmology: observations quasars: 

absorption lines - gravitational lensing - quasars: individual (Q1422+231) 



1. Introduction 

Much astronomical effort has been expended to investigate the structure of the universe 
on ever larger scales. Surveys of galaxies in emission, studies of intergalactic gas clouds 
in absorption, and searches for anisotropy in the cosmic microwave background indicate a 
universe approaching homogeneity on scales larger than a few hundred Mpc (e.g., Peebles 
1993). Velocity and density variations, smoothed over increasingly larger volumes, decrease 
steadily, in agreement with current cosmological n-body and hydro-simulations. 

In this paper, we investigate how the density and velocity fluctuations behave on the 
smallest scales. Even limiting ourselves to randomly chosen regions (which have densities far 
short of those encountered in stars and other highly collapsed objects) the answer is complex 
and must depend strongly on the local mean density. Large density fluctuations on very small 
(parsec) scales occur in the interstellar medium. The intergalactic medium (IGM), at least as 
observed at high redshifts, occupies an intermediate range between large scale homogeneity 
and small scale fragmentation. The combination of low density and photo-ionization heating 
by the meta-galactic UV background conspire to prevent catastrophic collapse and cooling. 
In the absence of discrete sources of momentum, heat, or ionization, any small scale structure 
below the Jeans mass should be smoothed out by the thermal pressure of the intergalactic 
gas. 

The nature of small scale structure in the IGM has received renewed interest with the 
advent of high resolution gasdynamic simulations of structure formation. It is possible to 
measure various cosmological and astrophysical quantities by comparing hydro-simulations 
of the general baryon distribution with observations of Lya forest absorption lines (e.g., Cen 
et al. 1994; Petitjean, Muecket, & Kates 1995; Zhang, Anninos & Norman 1995; Hernquist 
et al. 1996; Miralda-Escude et al. 1996; Wadsley & Bond 1996; Zhang et al. 1997; Croft et 
al 1997; Theuns et al. 1998; Bryan et al 1999; Dave et al 1999). In spite of the successes of 
the models in reproducing the observed Lya forest absorption, there is some concern that 
the finite mass resolution may limit quantitative conclusions as to the thermodynamics and 
ionization state of the gas. For example, the optical depth of an absorption line for a given 
value of the background density is approximately proportional to the dumpiness of the gas. 
Consequently, a locally denser gas could, to some degree, mimic a universe with a a larger 
baryon density and still satisfy the observational constraints (see Ranch et al 1997; Weinberg 
et al. 1997). Thus, estimates of the baryon density in the intergalactic medium depend on 
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either the simulations properly resolving those low- to intermediate column density structures 
in the IGM that contain most of the baryons, or on the absence of any significant small scale 
structure. Since the Jeans mass Mj varies proportional to p"^^"^ (oc (1 + z)~'^/^), there is 
always a redshift where the mimimum resolved mass exceeds Mj, so small scale structure 
present at high redshift could pose a problem. There are good reasons for concern that such 
structure may exist at some level. Velocity gradients of several kms~^, or appreciable density 
differences over distances of a few hundred parsecs in low column density gas, might imply 
that already at high redshift much of the matter is trapped in unresolved, small potential 
wells; alternatively, the gas may have been stirred by non-gravitational agents such as galactic 
winds and explosions (e.g., Couchman & Rees 1986). The widespread metal enrichment in 
the intergalactic medium shows that such galactic outflows must have occured. There have 
been several theoretical attempts trying to understand this result (e.g., Ferrara, Pettini, & 
Shchekinov 2000; Theuns, Mo & Schaye 2001, ; Cen & Bryan 2001; Madau, Ferrara, & Rees 
2001; Aguirre et al 2001), but little is known about the hydrodynamical consequences of 
these ancient winds for the IGM. Finally, it has been argued that even at high redshift many 
Lya forest absorption lines are formed in extended galactic halos (Chen et al 2000), in which 
case one should expect to see small scale hydro-dynamic disturbances in the gas caused by 
stellar feedback. 

The only way to test for the possibility of hidden small scale structure and of a universe 
grainier than envisaged by the simulations is to measure by how much the universal density 
and velocity fields fiuctuate on a given spatial scale. Differences in the absorption pattern of 
the Lyo; forest in closely spaced lines of sight to gravitationally lensed QSOs can be used to 
measure the dumpiness of the IGM on sub-kpc scales, about an order of magnitude smaller 
than can currently be resolved by cosmological simulations. Such observations have been 
used to estabhsh lower limits on the sizes of the Lya forest clouds (e.g.. Young et al., 1981; 
Weymann & Foltz 1983; Foltz et al. 1984; Smette et al. 1993,1995; Bechtold k Yee 1995; 
Petry et al 1998). Most of the above investigations have shown the difference between the 
lines of sight to be small. In the previous paper (Ranch, Sargent, & Barlow 2001, paper II), 
the second in a series of papers on the small scale properties of various gaseous environments 
at high redshift, we found that the metal enriched gas observed in CIV absorption systems, 
apparently produced in extended regions around redshift galaxies, shows signs of having 
been disturbed hydrodynamically on time-scales similar to those relevant for recurrent star- 
formation. This is consistent with the CIV gas clouds being the result of ancient galactic 
outflows. The CIV absorbing gas is mostly related to strong saturated Lya forest systems, 
and the question remained open as to how far the lower density IGM (causing the much more 
frequent, unsaturated Lya forest systems) has been affected. In the present paper, we use 
the Keck HIRES spectrograph (Vogt et al. 1994) to investigate the fabric of the Lya forest 
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proper on sub-kpc scales at high signal-to-noise ratio and with a velocity resolution sufficient 
to resolve the width of the Lyman lines and detect velocity gradients of a few hundred ms~^. 

The most suitable object known for this sort of study is the ultrabright, quadruple, 
gravitationally lensed QSO Q1422-I-231 {zem = 3.61; Patnaik et al 1992). This object has 
been observed for similar purposes previously at low resolution. Bechtold & Yee (1995) 
concluded from ground-based spectra that there are differences of typically about 50% in HI 
column density for beam separations on the order of 100 parsec. In contrast, Petry, Impey & 
Foltz (1998), from HST FOS spectra found very little difference both in column density and 
velocity between the various lines of sight. In section 2 we describe Keck HIRES (Vogt et al. 
1994) spectra of the A and C images of Q1422-I-231. The properties of this lens system are 
summarized on the CASTLES web site (http://cfa-www.harvard.edu/castles/B1422.html). 
The observations of images A and C of Q1422-I-231 and the methods used to reduce the 
data are discussed in section 2. Section 3 contains an account of the measurement of the 
differences in the Lyman a lines in the spectra of the two images. Our conclusions are 
summarized in section 4. 



2. Observations and data analysis 

Q1422-I-231 consists of four images. Images A (g= 16.92 magn.), B (16.77 magn.) and 
C (17.44 magn.) are almost along a line with C separated by 1.64" from A and 0.76" from B. 
Image D (20.56 magn.) is well separated from the rest. We chose to obtain HIRES spectra of 
A and C, using a 0.56" slit and only observing in conditions of excellent seeing. Much of the 
guiding was done manually and the position angle of the slit was chosen so as to keep image B 
out of the slit. In some exposures we guided on the outer edge of the image in order to ensure 
minimal contamination from B. The dataset is the same as described in papers I (Ranch, 
Sargent & Barlow 1999) and II (Ranch, Sargent & Barlow 2001). The spectra were reduced 
as described in Barlow & Sargent (1996). The individual A and C exposures were matched 
in flux to a template produced by adding up all A spectra, using polynomial fits typically of 
8th order per echelle order. Care was taken to include only points with fiux levels close to 
the continuum in the fitting regions, and to avoid absorption lines. This way, the large scale 
features in the spectra become identical (the information on large scale differences between 
the QSO continua is lost) but features on scales up to ~ 300-400 kms~^ are retained. Thus, 
differences among individual absorption line profiles up to that velocity scale survive our 
matching procedure and remain unaffected. The resulting spectra were normalized to a unit 
continuum using spline functions. Voigt profiles were fitted to the absorption lines using the 
fitting routine VPFIT (Carswell et al. 1991, see http://www.ast.cam.ac.uk/ rfc/ vpfit.html). 
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until the fit exceeded a minimum probability (1%) to produce by chance a x^ as high as 
the one attained. We departed here from previous such analyses in that the whole spectrum 
was fitted (in continuous pieces), without selecting significant absorption features by eye. 
We thus do not have to assign significance thresholds to individual lines, a procedure which 
is most meaningful for single unblended, strong absorption components, but which becomes 
uncertain in more complicated situations. Another problem, peculiar to relatively high 
signal-to-noise ratio spectra such as those obtained with HIRES, is that it is hard to obtain 
satisfactory fits merely by superimposing Voigt profiles, even if Voigt profiles are a good 
physical model. When S/N ratios exceed a few tens, the errors in the fit are invariably 
dominated by systematic uncertainties in the placement of the continuum level. Therefore, 
it was decided to treat the continuum in every fitting region (average length 28 A) as a free 
parameter. This resulted in a mean drop of the continuum level over the region [4747, 5630] 
A (i.e., between Ly/5 and Lya emission, where the Lya analysis in the present paper was 
performed) to 99.3 % of the original level, with a standard deviation of 2.7% of the original 
continuum level, i.e., the total continuum level remained virtually unaffected by introducing 
locally free continua, but the local continua were adjusted typically by 2.7% up or down, to 
obtain an optimal fit. 

Since the A and C images are so close to B we have to be concerned that blending of the 
spectrum of an image with light from the other images might have reduced the differences 
seen between the spectra. However, our spectra show strong column density differences 
between the A and C images in some low ionization metal lines in the Lyman a forest 
(see paper I) which would have been significantly reduced if spillover of light from different 
images had been a serious problem^. Our technique of matching the two spectra to a common 
template as described earlier does take out the large scale fiuctuations but not those among 
individual lines. Again, the differences seen in metal line systems suggest that this is not a 
problem. Finally, while they may affect the measured column densities, none of these effects 
could conspire to reduce any actual intrinsic velocity differences between the lines of sight to 
values as tiny as those found here. 

The lowest redshift covered by the selected spectral region, z = 2.905, corresponds to a 
transverse separation Ar ~ 270 pc. Throughout the paper beam separations are computed 
for Hq = 50 kms~^Mpc^^ and go = 0-5. The redshift of the lens is taken to be z = 0.338 



^we have also obtained as yet unpublished spectra of the Icnscd QSO UM673 with HIRES. In this case 
the images A and B have a separation of 2.2" so there is little risk of mixing the light from different images. 
The spectra have significantly lower signal-to-noise ratio than the data on Q1422+231 discussed here, but 
the absence of differences between the low column density Lya forest lines in UM673A,B is consistent within 
the errors with the very small differences found in the present paper. 



-6- 



(Kundic et al. 1997). 



3. Differences between the lines of sight in the Lya forest region 

Here we describe three main approaches to measuring differences between the Lya 
forests in the two adjacent lines of sight: a global cross-correlation analysis; a search for 
column density and velocity differences among individual absorption lines; and a pixel by 
pixel measurement of relative fluctuations in the optical depth between the lines of sight. 



3.1. Global correlation analysis 

We can study global differences in the Lya forest region by measuring the cross- 
correlation function, ^cc, over the total, useable length of both spectra. We define this 
quantity by 



icci^v, Ar) 



< {F^{V)- <Fr>)- (Fr+Ariv + At;)- < F,+Ar >) > 
'< (Friv)- < Fr >)2 > • < (Fr+Ar(t^ + Av)- < F^+Ar >)' > 



Here F^ and -Fr+Ar are the pixel flux values of the two spectra, separated by Ar on the 
plane of the sky. The velocity coordinate along the line of sight is v (where dv = cdX/X), 
and Av is the velocity lag. The averages are taken over the whole velocity extent of the 
spectrum. For Ar = we get the usual auto-correlation function C,cc{Av, 0), while for Av = 
we have the cross- correlation as a function of transverse separation only. The function is 
defined so as to satisfy .^cc(0,0) = 1. With large scale velocity correlations (> 1000 kms^^) 
expected to be absent or weak (Sargent et al. 1980), the autocorrelation function (on scales 
~ 100 kms~^) mostly measures the Lyct line width and the weak small scale clustering of 
Lya forest systems (e.g., Webb 1987, Ranch et al. 1992). We apply the correlation analysis 
to most of the spectral region between Lyman f3 and Lyman a emission, [4737, 5630] A. 
Only regions [4747, 4875], [4895,4997], [5003,5221], [5227,5409], [5419,5473], and [5480,5630] 
A were used in the analysis in order to avoid contamination by known metal systems. These 
regions are shown in fig. 1; the hatched parts were excluded. The resulting mean redshift z = 
3.26036 of the remaining sample corresponds to a mean beam separation Ar (=119 h^Q^pc). 
A small portion of the spectrum is shown in fig. 2. The spectrum of image C (dotted line) 
is overplotted on that of A (solid line), to illustrate how closely the spectra resemble each 
other. 
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The function ,^cc(Af , Ar) is shown in fig. 3. In particular, we obtain the "zero-lag" 
cross-correlation function: 

^{Av = 0; Af = 119 h^^ pc) = 99.48 %. (2) 

Thus, over a mean transverse distance of 120 pc the appearance of the Lya forest changes 
very little, and even the small difference observed may be due to noise or systematic errors. 



3.2. Comparison between individual absorption lines 

We next investigate how properties of individual absorption systems change with the 
separation on the plane of the sky. We can quantify differences Ax = Xi — X2 in any physical 
property x of absorption systems common to both lines of sight 1 and 2 by means of the 
distribution function of Ax. (For example, consider the difference between the redshifts of 
the lines, Zi — 22). In an isotropic universe, the mean. Ax, vanishes, so the information 
is contained in the shape of the distribution of differences, as characterised, for example, 
by the observed scatter of the distribution, o-(Ax). The width of the intrinsic distribution, 
aint{Ax), differs from the observed one due to the combined measurement errors, o'{xi) and 
0"(x2), and any systematic error CTgys resulting from the occasional erroneous pairing of two 
unrelated absorption components. The quantity o"<(Aa:) = y^a{Axy — a(xiy — o'(x2)'^ is 
then a strict upper limit on the intrinsic scatter, aint{Ax). 

Here we concentrate on the distribution of velocity and column density differences, Av 
and A log A^. Inspection by eye shows that there are only minute difference between the 
spectra of the A and C images. Accordingly, it is legitimate to associate each Lyman a 
forest line in the spectrum of component A with a corresponding line in C Thus, after 
profile-fitting the spectra, we suppose that an absorption component in spectrum C belongs 
to the same underlying physical structure as a line in A if its redshift and column density 
fall within a given redshift and column density window centered on a similar component in 
the A spectrum. 

We make the assignments totally automatically in order to preclude human prejudice. In 
order to avoid too many random associations we restrict the allowed range of column density 
differences when measuring the relative velocity distribution, and restrict the allowed range 
of velocity differences when comparing the column densities. Specifically, the redshifts had to 
agree to within ± 30 kms~^ (there is no real prejudice as this is much larger than the actual 
scatter), while column density differences jlogN^ - logNcI < 1 were deemed acceptable. 
In order to avoid including ill-constrained measurements and spurious cross-identifications, 
all lines included were required to have relative Doppler parameter errors, a{b)/b <3 and 



column density errors a{logN)/logN<3. Doppler parameter values had to fall in the interval 
[5, 300] kms~^ to exclude spuriously low values due to noise spikes and high values due to 
unresolved blends. All error estimates given below are 1-sigma deviations derived from the 
output of the line profile fitting program. Finally, to get rid of obvious mis-associations, one 
iteration of a 4-0" clipping algorithm was performed. Column density differences larger than 
4 standard deviations of the mean of the distribution of column density differences where 
rejected. The analysis was limited to column densities 12< logA^ < 14.13, the lower bound 
being close to the detection limit, whereas the upper limit (14.13) corresponds to a residual 
flux of about 2.4 % (for a Doppler parameter of 28 kms^^, typical of an unsaturated 2; ~ 3 
Lya forest line), a value just at the top of the linear part of the curve of growth, but far 
enough away from saturation to permit us to see significant differences between the central 
optical depths. This column density regime is most relevant for the study of gravitational 
structure formation since the selected absorption lines probe gas overdensities still close to 
the linear regime of gravitational instability (Ap/p < 15) on the scale of individual galaxies 
(e.g., Miralda-Escude et al. 1996). 



3.2.1. velocity differences 

Fig. 4 shows the observed histogram of velocity differences, Av = va—vq^. The top plot 
of fig. 5 shows the observed standard deviations of the distribution of velocity differences 
between the two lines of sight (solid dots), versus the largest measurement error in redshift 
for lines remaining in the sample. The open diamonds show the predicted width of the 
distribution of velocity differences expected from statistical measurement errors alone. It 
was computed as the quadratic sum of the 1 — a error estimates produced by the Voigt 
profile fitting routine. To search for intrinsic scatter we discard those velocity differences 
with large errors. As the most uncertain points are eliminated, the observed width of the 
distribution should decrease, and any intrinsic scatter should emerge. We see that, the 
further one goes to the left along the abscissa, the line pairs remaining in the sample have 



®The original distribution did not have a zero mean but showed an mean shift of spectrum C by 0.96 
kms^^. This shift amounts to sUghtly less than half a pixel or a quarter of the slit width, fts size and 
the circumstances of the observations make it most likely that it was caused by systematic variations in 
the position of the target on the slit: the targets were only observed in exceptional seeing to minimize 
contamination of the desired image by the light from other images. Consequently, the seeing may often have 
been better than the 0.574" slit width, leading to a dependence of the wavelength zero point on the actual 
position of the target in the slit. This uncertainty may have been turned into a systematic shift by our 
attempts to avoid spillover. For the analysis below and for all the plots from Fig. 4 on, the above shift was 
subtracted from the velocities of the C spectrum. 
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smaller and smaller maximum (and mean) redsliift errors and the more well defined are 
the measurements. In other words, if the bulk of the scatter is caused by measurement 
error, we should expect to see the observed total scatter a{Av) decreasing in step with, and 
closely tracking, the measurement errors as increasingly stringent error limits are applied 
and more and more data points are discarded from the distribution. Aside from a certain 
scatter (probably caused by mis-matched components) this behaviour is indeed observed in 
fig.5. (The discrepancy at large values of the abscissa between the width of the predicted 
distribution based on measurement errors alone, and the observed distribution, is an artifact 
of the range of parameters we imposed earlier. It is caused by our admitting only velocity 
differences < 30 kms~^ into the sample. This is perfectly justified given the range of the 
large majority of values in figure 4, but it does get rid of some large outliers which result 
from automatic mis-identifications of line pairs.) 

To demonstrate the sensitivity of this method, we plot (for illustrative purposes) the 
observed a{Av) a second time (dashed line in the top panel of fig. 5), but now with an 
additional random, Gaussian velocity jitter of a=1.5 kms^^ imposed on each line pair, i.e.. 



the distribution of pair velocity differences has now a width of \/a'^{Av) + (1.5 kms~^y. On 
the left half of the diagram this jitter clearly dominates the observed width, demonstrating 
that any intrinsic velocity differences of that size would have been seen easily and that any 
actual differences must be considerably less than 1.5 kms~^. 

There are 259 line pairs with individual redshift errors o'{z)< 4xl0~^ with an observed 
total cr{Av) =1.7 kms~^, and still 83 pairs in the bin with the smallest errors, (j{z)< 1 x 10"^ 
(o'(Af) =0.6 kms~^). The upper limits o"<(Af) on the intrinsic aint{Av) (where the mea- 
surement errors have been subtracted in quadrature from the observed width) are formally 
0.4 and 0.2 kms~^, respectively. Given the scatter between the measurement errors and total 
a{Av) visible in the plot, the intrinsic distribution is probably consistent with zero width. 
We can compare these values with the velocity differences expected from the Hubble expan- 
sion alone. If Lya clouds expand with the Hubble flow the projected velocity differences 
between two lines of sight, (modelling the geometry by an expanding slab intersecting both 
lines of sight at an angle a) , would be on the order of Avh/ tan a = H{< z >) r ~ 50 ms~^, 
where H{< z >= 3.26) = 440/i5o kms~^Mpc~^ was adopted as the Hubble constant at the 
mean observed redshift. This is well below our detection threshold and our null result is 
consistent with the clouds following the Hubble flow. 

It is tempting to try using the measurement of (x{Av) to get a rough idea of the energy 
deposited in the gas in the form of turbulence. In paper II we had applied the Kolmogorov 
relation between the structure function B{r), and the energy input rate e at scale r, 

B{r) ^ [vi{r') - V2ir")]^ ^ (e r)^/^ (3) 
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(e.g., Kaplan & Pikelner 1970; here f i(r') and V2{r") are the velocities along the line of sight, 
and r is the spatial separation between the lines of sight, and the average is taken over all 
absorption line pairs with mean transverse separation r) to get an approximate estimate of 
the energy input rate e in CIV absorbing clouds. These gas clouds correspond to strong, 
usually saturated Lya forest systems and overdensities higher than a few tens. This approach 
could be justified then as we had indeed observed turbulent energy on a wide range of scales 
and had found a rise of B{r) with r; moreover the relevant time scales for energy dissipation 
were short, so turbulent energy must have been propagating down over a range of scales 
close to the epoch of observation. Now we are applying this scheme to the lower density 
gas giving rise to unsaturated low column density Lya lines, which has densities within a 
factor of 10 - 15 of the mean density of the universe, and, at average, must be situated much 
further away from sources of stellar energy. Here it is more likely that any energy injection 
by stellar feedback is intermittent or happened only once long ago so a steady state (one of 
the conditions for the validity of the Kolmogorov law) may not exist (or no longer exists). 

With this warning in mind we proceed to calculate e for our present low column density 
Lya forest sample. The quantity e is still a measure of the energy fiow on a given scale, 
but is not necessarily related to the current primary energy input on larger scales (e.g., by 
explosions or winds). If we identify our measured variance a'^{Av) with B{r), we can solve 
for the energy transfer rate^ at scale r: 



e ^ -. (4) 

r 

Adopting the above value (Ti„((Af ) = 0.4 kms^"*^, and f = llOh^Q pc, 

e ~ 2 X 10~^/i5o erg g^^s"^ (5) 

Thus the intergalactic medium giving rise to the unsaturated parts (i.e., most) of the Lya 
forest at z=3.26 is experiencing energy transfer at a rate three to four orders of magnitude 
less than the gas giving rise to the CIV absorption systems dealt with in paper 11^, and six 
to seven orders of magnitude less than the values in the Orion nebula, for example. 



^Because of the uncertain geometry of the clouds and our taking the known projected transverse separation 
instead of the (unknown) actual 3-diinensional separation between points (see the discussion in paper II) 
this approach can only provide a rough, order of magnitude estimate of e. 

^but note that the CIV turbulence was measured on a three times larger scale, 300 pc, so the discrepancy 
may be less. 
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3.2.2. column density differences 

A scatter plot of the column densities in A and C is given in fig. 6. The scatter is 
larger in the C image because of the lower signal-to-noise ratio there. Some of the outliers 
may come from the ambiguity of the Voigt profile models (which occasionally may lead to 
similar absorption complexes being modelled automatically by two rather different sets of 
Voigt profile components), while others may be genuine differences due to unidentified metal 
lines. The distribution of column density differences is shown in fig. 7. Two components 
were considered to correspond to each other, when they were nearest neighbors, lay within 
30 kms~^ of each other, and had column densities differing by less than 1.0 dex. Fig. 8 
shows the scatter versus measurement error diagram for column densities, similar to Figure 
5. For the best defined samples (those with measurement errors around 5xl0~^ in logiV) 
the corrected upper limit to the intrinsic scatter is o"<(AlogiV)~ 2 x 10~^. We note that it 
makes sense to discuss an 'optimum sample' because discarding those data points with the 
largest measurement error leaves better measurements but also reduces the sample size and 
thus the statistical significance. Initially throwing away the points with the worst error will 
improve the precision with which sample mean and variance are known, but the precision 
of these statistical values decreases with the square root of the number of data points, yM 
and if enough points have been lost the penalty gets larger than the gain. In the appendix 
we have defined a function of merit, E, which is positive as long as it is beneficial to discard 
data points, and crosses zero towards negative values if the sample becomes to small. The 
bottom panels of figures 5 and 8 show this function. For example, from the column density 
case (fig. 8) we find that the sample has been cut down to optimal size by retaining only 
data points with a maximum column density error < 0.05. 

Most interesting from an astrophysical point of view is whether we can put limits on 
the total baryon density fluctuations on the scale of the transverse separations of the two 
lines of sight. The relation between the HI column density N(HI) and the baryon density for 
Jeans-smoothed IGM regions can be written as 

N{HI) oc pta, 

with a = 1.37 - 1.5 (Schaye 2001). For this power law the distribution of column density 
differences AlogA^ (= logA^^ ^ ^og Nb) can be transformed by simple scaling with 1/a = 
0.73 (we use the smaller value a = 1.37 to establish the most conservative upper limit) 
into a distribution of density differences, and the variance of the logarithmic baryon density 
becomes 

{{A\ogpf)<a-m\ogNA-\ogNBf) (6) 

Adopting (conservatively) the pair sample with mean column density measurement errors 
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below 0.05 dex (153 lines; see fig. 8) we arrive at a value of 2.5 x 10^^ for the corrected 
"intrinsic" scatter and 4.3 x 10^^ for the uncorrected, observed scatter in AlogA^. As the 
figures show, other choices of sample size give similar values. Then we get 



(Alogp) ) < 1.8 X 10"^(or 3.1 x 10"^ uncorrected), (7) 

for the typical logarithmic change in density on a scale of 0.110 h^Q kpc, i.e., the RMS 
fluctuation in the baryon density is less than about 3.1 percent. From the close agreement 
between the width of the observed distribution and the expected width from measurement 
errors alone, it is clear that the fluctuations are dominated by the measurement errors and 
that the intrinsic fluctuations are consistent with being zero at our level of precision. 



3.3. The filling factor of 'disturbed' regions in the universe 

Here we address the question as to how much of the universe (by volume) is hydrody- 
namically disturbed, for example by supernova explosions or galactic winds along the two 
lines of sight. We consider a region in space at a given redshift z as 'disturbed', if the optical 
depth r for HI Lya absorption differs by more than a certain amount between the two lines 
of sight. 

This is easy to measure unless the lines are saturated. Many of the saturated systems 
are associated with multi-component metal absorption systems. These systems exhibit sub- 
stantial variations in the absorption pattern across multiple, close lines of sight (e.g. Lopez et 
al 1998; paper II), but the component structure and any differences in the column densities 
are almost always completely obscured in the strong, saturated HI Lya absorption blends. 
Thus, to get a conservative upper limit on the fraction of disturbed regions we assume that 
saturated regions in the Lya forest are always disturbed, i.e., we treat them as if they had 
infinite optical depth differences. 

Fig. 9 shows the distribution of the difference in optical depth At =\ ta — tq \ between 
the two lines of sight A and C. The four different histograms show the results for various 
column density ranges^ (the column density in spectrum A was used to determine whether a 
pixel pair fell into a certain range): the solid line is for the whole spectrum, the dashed one 
for weak lines with 12.0 < logA^ < 13.0, the dash-dotted one for 13.0 < logA^ < 14.0, and 
the dotted one for stronger, but still unsaturated, lines. (A column density logA^ = 14.13 



®the column densities quoted here are calculated by treating the observed optical depths as if they were 
central optical depths of Gaussian lines with Doppler parameter 28 kms~^. 
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corresponds to a residual flux of 2.4% for the center of an absorption line with a typical 
Doppler parameter 6 = 28 kms~^). In the first histogram (solid line), saturated pixels 
(for practical purposes those with logA^ > 14.13), did not provide any information and 
were assumed to have infinite optical depth differences Ar; thus they do not appear in this 
histogram in the lower bins. The relative proportions of the spectra in those four column 
ranges was 100%, 36.6%, 36.4%, and 17.4%, respectively. To reduce contributions to the 
optical depth differences from noise, both spectra were smoothed with a box filter on the 
scale of typical line broadening (28 kms^^). 

We found that 76.8%, 98.7%, 74.7%, and 46.7%, respectively, of all the pixels in the 
four column bins had their optical depths differ between the two lines of sight by less than 
Ar = 0.05 (i.e., they contribute only to the first bin of fig. 9). It is remarkable that of the 
36.6% of all pixels which are in the column density range 12.0 < log A^ < 13.0, 98.7% agree 
to within r < 0.05. Apparently, the low column density Lya forest is virtually undisturbed 
by any effects capable of producing differences in density over a couple of hundred parsecs. 
The differences in optical depth increase with increasing column density as one would expect 
if higher column density gas is more closely associated with galaxies, which are the likely 
origin of any hydrodynamic disturbances. 

The overall fraction of the Lya forest with Ar < 0.05 irrespective of column density is 
1 — IloS = 76.8%. Since we have counted regions close to or beyond saturation as having 
infinite Ar (see above), JloS is a strict upper limit on the line-of-sight filling factor of the 
Lya forest with absorption from processed gas capable of producing a disturbance of that 
magnitude. 



3.3.1. an upper limit to the filling factor of galactic wind bubbles 

The upper limit fios ^ 0.23 on the fraction of disturbed pixels in the Lya forest can 
be used in a crude way to estimate the number density and volume filling factor of galactic 
outflows from their effect on the density field of the IGM, in sofar as such variations appear 
as differences between the HI column densities (Lya forest opacities) between two close lines 
of sight. The principle will be illustrated here for a model of windbubbles expanding into 
the intergalactic medium (fig. 10). 

We assume that dense shells of swept up material occasionally cross the double lines 
of sight to a lensed background QSO. It can be shown (see Appendix B) that even under 
unfavorable conditions the shells should have HI column densities high enough to be seen in 
absorption in a high signal to noise ratio QSO spectrum like the one of Q1422+231 discussed 
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here (N(HI)> lO^^cm"^). The passage of the shell is very likely to disturb the resemblance 
between the absorption lines in the two lines of sight (Fig. 10). Ideally, the shell should 
produce two absorption systems at its two points of intersection with the lines of sight. We 
assume that there will be measurable differences in the absorption lines between the lines of 
sight as long as new shell material keeps crossing the lines of sight faster than the differences 
can be smoothed out by pressure waves travelling at the sound speed. Or in other words, a 
shell is no longer detectable for when the speed of sound c exceeds the transverse velocity of 
the shell across the lines of sight: 

^b < c, (8) 

where b is the impact parameter of the lines of sight relative to the center of the shell, and 
R and R are the radius and expanding velocity of the shell. 

To convert the line-of-sight filling factor into a volume filling factor a particular model for 
the disturbance and its observational signature is needed. For illustrative purposes we adopt 
the supershell model of MacLow & McCray (1988), i.e., we assume that the disturbances 
are wind-blown shells escaping from high redshift galaxies. For a given line-of-sight filling 
factor, the radius R, velocity of the shell R, volume filling factor f^, and the space density 
Hgai of galaxies with a wind escaping into the IGM can be written in terms of the mechanical 
luminosity L38 (in units of 10^^ ergs/s), the density of the ambient medium n_5 (in units of 
10^^ cm^^), and the age of the wind bubble ty (in units of 10^ years). These relations are 
given in Appendix B. 

Figures 11 and 12 show the dependence of the upper limits on Ugai and /^ for winds 
of various 'strengths' (i.e., with the ratio L^^/n^^ as the distinguishing parameter) corre- 
sponding to the upper limit JloS < 0.23 on the fraction of disturbed pixels in the Lya 
forest. 

The model is highly oversimplified in that it assumes spherical bubbles with a single 
age and a single combination of L38 and ngau propagating through a homogeneous medium 
(and the expansion of the universe is not taken into account), but it contains several features 
likely to be relevant to more realistic future models (figs. 11 and 12): 

1. The number of galaxies with winds originating close to the epoch of observation {z = 3.26) 
is essentially unconstrained by an observation of the line of sight filling factor (fig. 11). This 
is because the shell radii are so small that even a large number of galaxies with winds would 
produce a negligible total absorption cross section. 

2. Weak winds (low luminosity and / or large density of ambient medium, i.e., L^s/n-5 < 1), 
because of their limited lifetime, escape detection if they originated early (e.g., before zsf ~ 4 
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for a L38/n_5 ~ 0.1 shell, or before Zsf ~ 7 for a L^s/n^^ ~ 1 shell; fig. 11). 

3. Strong winds (Lsg/n-s > 10) remain detectable even if they started a Hubble time ago 
(fig. 11), and the upper limit on their numbers becomes stronger with increasing redshift 

^SF- 

4. The volume filling factor (fig. 12) is most tightly constrained for recent rather than old 
winds, independent of strength. 

5. Our observation cannot exclude that the universe is filled with ancient, weak wind bubbles 
because they would no longer be detectable with the present method by redshift z ~ 3 
(see item 2 above). However, since these weak winds do not reach far out from galaxies 
and galaxies are not distributed randomly it is very unlikely that such bubbles could fill a 
substantial fraction of the voids. Strong winds (Lsg/n.s > 10) are much better constrained 
because of their long lifetimes: if winds of strength L^s/n^^ ~ 1000 started as early as 
Zsf ~ 10, they could be filling up to about 40% of the volume of the universe by redshift z 
= 3.26. If such strong winds arose at Zgf ~ 4 they may fill up to 18% . 



4. Conclusions 

We have applied various methods to search for and measure density and velocity struc- 
ture on 100 pc scales in the intergalactic medium aX < z >= 3.26. The smallness of the 
differences in column density and velocity between the lines of sight confirm and extend the 
results of Smette et al (1992,1995), and Petry, Impey & Foltz (1998). In the unsaturated 
regions of the Lya forest RMS differences between the projected line of sight velocities over 
120 pc do not exceed a few hundred meters per second. The RMS total density differences 
are inferred to be less that 3%, and are basically below the detection limit. Visible differ- 
ences between the lines of sight do occur occasionally but are localized and seem to be mostly 
limited to high column density gas. Strong differences are known to occur from studies of 
the higher column density metal absorption lines (papers I and II), but these differences 
rarely show up in the Lya region proper because of saturation. Our results appear broadly 
consistent with predictions by Theuns, Mo & Schaye (2001) who model the impact of stellar 
feedback on the appearance of the Lya forest and find that stellar feedback will mostly affect 
the higher column density, Lyman Limit absorption systems. The suggestion (Chen et al 
2000) that most absorption lines with a rest frame equivalent width > 0.32A (corresponding 
to N(HI)~ 6 X 10^^ cm ~^ for a Doppler parameter b = 28kms~^) are related to galactic halos 
(where we might expect to see hydrodynamic disturbances) is not contradicted by our results 
because close to 90% of the Lya forest lines above our chosen column density threshold (10^^ 
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cm ^) have lower column densities. 

Our results are consistent with unsaturated Lya systems being mostly featureless at the 
few percent level on scales smaller by an order of magnitude than state-of-the art cosmological 
hydro-simulations can currently resolve. Apparently, the finite resolution of the simulations 
does not lead to a significant underestimate of the dumpiness of the baryon distribution, 
for densities typical of the unsaturated Lya forest, i.e., overdensities up to a factor 10-15. 
However, our observations only provide a snapshot of the properties of the Lya forest at 
z ~ 3.3 and cannot exclude the existence of earlier non-gravitational processes in the IGM 
if, by the time of observation, their hydrodynamical traces have been obliterated (see also 
below). 

Notwithstanding doubts about the applicability of the Kolmogorov law to the low den- 
sity IGM, we have tentatively derived the rate of turbulent energy input into the low 
density IGM. The turbulent energy in the unsaturated low column density Lya clouds 
{log N{HI) < 14.13) appears to be lower by 3-4 orders of magnitude than that in the 
stronger, CIV absorbing clouds dealt with in paper II. Apparently, the CIV gas has been 
affected more recently (and thus more frequently, given it is unlikely that there is a sharply 
defined "epoch" of star formation) by stellar feedback than the general Ly« forest. 

The fraction of the spectra that differ by more than 5% in optical depth with respect 
to the other line of sight, considering all saturated regions as being disturbed, was measured 
for various column density ranges. The results imply a limit on the relative importance of 
hydrodynamic disturbances in the presence of the restoring effects of gas pressure. From the 
upper limit to the line of sight filling factor {fios < 23%) we find upper limits to the volume 
filling factor and space density of wind-producing galaxies as a function of the mechanical 
wind luminosity, the density of the ambient IGM, and the starting redshift of the wind. The 
presence of short lived, weak winds (low luminosity or high ambient density) is only poorly 
constrained. The other extreme, long-lived superwinds with Lag/n-s ~ 1000 arising as early 
as redshift 10, could fill up to 40% of the volume of the universe by redshift 3.3, but not 
more than 18% if starting at redshift 4. Currently there is no observational evidence for the 
existence of 2; ~ 10 superwinds, but closer to the epoch of our observations strong outfiows 
have been seen in high redshift galaxies (Franx et al. 1997; Pettini et al 2001), so there is 
the possibility that a significant fraction of the volume of the universe may have been stirred 
by winds; our results appear to show that it is not a dominant one. 

MR is grateful to NASA for supporting the initial phases of this work under HF- 
01075.01-94A, and to the NSF for grant AST-0098492. The work of WLWS was supported 
by NSF grant AST-9900733. We thank Martin Haehnelt, Joop Schaye and an anonymous 
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5. Appendix A 

Our criterion for the optimum choice of the maximum permissible measurement error 
(chapter 3.2) is as follows: we would like to omit m of the n measurements such that the 
mean measurement error decreases from o"crr {n) to ae„ {n — m), as long as this reduced error 
outweighs the decrease in the number of data points available, and continues to lead to 
a reduction in the error of the determination of the total width of the distribution to be 
measured, Ototip)! \f^ -^ '^totiP' — m)/ \/n — m; in other words, we throw away the data 
points with the largest errors until the uncertainty in the width of the distribution does not 
improve anymore. Data points should be continued to be omitted as long as 

E = {n-1 + S{m)) ^Un)- ^^^^^^ ^^^.(n - m) - a^H S{m) > 0, (9) 

n — m 



where 



^/ N v^ "nxn — 1 , ^ 



For m = 1, we obtain the simple relation 



^tot(«) 



a2,,(n) - a2^,(n - 1) > -^Hi^ (11^ 

n 



The function E is shown for the distributions of the velocity and column density dif- 
ferences in the bottom panels of figs. 5 and 8. The optimum value of E occurs where E 
changes from positive to negative values. 



6. Appendix B 

This appendix addresses two questions: (1) would a superb ubble or wind blowing into 
the IGM at high redshift produce detectable HI absorption, and (2) can we produce limits 
on the volume filling factor of extragalactic wind bubbles from the fraction of the Lya forest 
spectrum disturbed by such events ? 



6.1. Detection of HI absorption from dense shells in the IGM 

We assume that the dynamics of the wind shell can be described by the superbubble 
model of Mac Low & McCray (1988). Then the dependence of the radius R on the mechanical 
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luminosity L = L^s x 10^^ erg s ^, the time t = t^ x 10^ years, and the number density of 
protons n = n_5l0~^ cm~^ in the ambient IGM, is given by 



R = ly^nZllhy^ 2670 pc, (12) 

and the expansion velocity is 

R = L^^^-J^hj^'^ 157 kms"\ (13) 

The total hydrogen column density seen radially outward through the shell is 

Ns{H) = UsAR = — = 2.75 X lO^^L^f n^'t^^'cm-^ (14) 

The gas density Ug in the shell is given by (Weaver et al 1977): 

ns = n — , (15) 

c 

where c is the sound speed in the undisturbed IGM and Cg is the sound speed in the shell. 
The thickness of the shell, Ai?, is 

^^^iV,^^^ E.^ - 14.4 ./^:4;r<'° PC. (16) 

The sound speeds are measured in units of 20 kms^^. The density in the shell is then 



r4r2/5„3/5 4/5 2 
Us = 2 '^^ ■ y^') 



It is not clear how exactly the gas is ionized. If it were in photoionization equilibrium with 
an ionizing UV background radiation field with intensity J = 10~^^ erg cm~^s~^IIz~^sr~^ 
the neutral fraction would be approximately (Bergeron & Stasinska 1986): 

^(^^) - 1.67 X 10-^ ^ "^ 



n{HII) Vl0-2cm-3. 

Thus the HI column density of the shell (from equations 14 and 18) is 



n{HI) 



.12 ^7/5^3/5,7/5-2 



^2 ^2/5.-2/5,-4/5 



^'^^^^ = ^imif'^^^ = ^-^ "" ^° ^-5 ^3ri7 c-2'0 [1 + 0.0162c>1^5U3t^%-^^°J cm-^(19) 

which is well detectable in absorption in our data. The actual absorption is likely to be 
stronger (the bubble certainly sweeps up material with higher than average density (n_5 ~ 1 
at 2; ~ 3) in the vicinity of the galaxy from which it originates. 
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6.2. Expanding shells as revealed by small scale differences in the absorption 

pattern 

A dense shell crossing two adjacent lines of sight separated by a distance Ar will produce 
differences in the column densities as long as the time it takes the shell to traverse the lines 
of sight is shorter than the sound crossing time Ar/c (= 5.8 x 10^ yr, in the present case). 
From simple geometric arguments, 

^b > c, (20) 

where b is the impact parameter of the lines of sight relative to the center of the shell. Thus 
the shell is detectable for b > bmm = cR/R, and the cross-section a of the shell for producing 
absorption lines with differences between the lines of sight that have not been ironed out by 
pressure waves is 

,2 ' - <- 



a = 7ri?M 1 - — 1 , for i? > c (21) 

0, for R<c. (22) 

The 'footprint' of the shell absorption in the Lya forest spectrum is expected to consist of a 
pair of absorption lines, and we assume here that each of these is broadened with a velocity 
width Avfwhm = 46 kms^^ (Doppler parameter 28 kms^^, which is typical for a photoionized 
Lya line in the z=3 IGM). Conservatively, we count only the sum of the widths of the two 
lines as contribution to the disturbance in the Lya forest spectrum as it is not clear whether 
the interior of the bubble would produce detectable structure between the lines of sight (it 
is most likely to produce a clearing in the forest without small scale features). 

When expressed in spatial coordinates (at z=3.26) the combined velocity widths of 
the two lines corresponds to a proper distance AL = 2AvfwhmH{z)~^ ~ 206/;,j^o^kpc. The 
fraction of the Lya forest spectrum disturbed by a coeval population of galaxies with space 
density Ugai, producing identical, spherically symmetric, expanding wind shells with radii R, 
and expansion velocities R would be 

has = ngai^AL = 2n,ainR'(l-^)^^^^, for R > c, (23) 

whereas the volume filling factor for the wind bubbles is simply 

f.=n,ai^R' = lhosR(l-t) '-^^^. (24) 

3 3 \ R^J Avf^hm 

Using the Mac Low & McCray superbubble model to get a first crude estimate the quantities 
ngai and /^ can be expressed as functions of the line of sight covering factor fios^ the time 
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ty or redshift Zsf when the wind (or star forming episode) started, and the ratio between the 
mechanical luminosity and the gas density in the IGM, Lsg/n-s. The upper limits on Ugai and 
/„ corresponding to the upper limit f^oS ^ 0.23 are shown as a function of Zgf for various 
values of the parameter Lsg/n.s in figures 11 and 12, respectively. The transformation 
between time and redshift assumes Qm = 0.3 and Hq = 50 km s~^ Mpc~^. 
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Fig. 1. — Normalized Lya forest spectra of the QSO 1422+231 A image in the wavelength region 
between Ly/3 and Lya emission. The shaded areas were excluded from the cross-correlation analysis 
as not belonging to the forest or as being contaminated by known metal absorption lines. The 
useable Lya forest region extends over the wavelength range [4747, 5630] A or redshift range 
z = [2.90,3.63]. The maximum beam separation at the short wavelength end of the spectrum is 
270 /igQ pc. The slowly varying function at the bottom is the 1 — a error. 
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Fig. 2. — Enlargement of part of the spectrum of Q1422+231. The spectra of the A (sohd hne) and 
C (dotted hne) images are plotted on top of each other. There is hardly any difference between the 
spectra, with the exception of a narrow line (which is largely absent in the spectrum of the C image) 
and some lower column density fluctuations near 5226 A. This line and the sharp stronger one at 
5225 A can be identified with a SilV 1394A interloper from an absorption system at z = 2.74889. 
The other absorption features are plausibly attributed to HI Lya. The slowly varying function 
plotted at the bottom is the 1 — a error of image A. 
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Fig. 3. — cross-correlation function r(< d >,Av), for a mean beam separation < d > 

pc, between the Lya forests in the Q1422+231 A and C image, for the spectral ranges shown in 

figure 1. 
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veloclty difference Av [km/s] 

Fig. 4. — Distribution of the velocity differences Av (= va - vc) for pairs of absorption components. 
Two components were considered to correspond to each other, when they were nearest neighbors, 
lay within 30kms~^ of each other, and had column densities differing by less than 50 %. 
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Fig. 5. — Top: RMS difference between the velocities of individually fitted absorption components 
in the two lines of sight, versus the maximum measurement error permitted for the redshifts in 
the sample used. Solid dots give the observed width of the Av distribution, open diamonds the 
width expected on the basis of statistical errors in the velocity determination alone. The dashed 
shows the width of the observed distribution, convolved with a Gaussian scatter with a = 1.5 
kms~^. Bottom: function E giving the merit of increasingly getting rid off the data points with 
the largest error (s.a. fig. 8). The value of E shows the trade-off between omitting measurement 
points with errors above a certain value, and the consequent reduction in the size of the sample, 
versus the maximum measurement error accepted in the sample. Positive values of E imply that 
it is beneficial to continue throwing away points with the largest errors, negative values imply a 
loss of information. The zero crossing corresponds to the optimum choice. The meaning of E is 
explained in the text and the Appendix. 
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Scatter diagram of the logarithmic column densities in the A and C spectrum. 
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Fig. 7. — Distribution of the logarithmic column differences A log N (= log Na - log Nc) for pairs 
of absorption components. Two components were considered to correspond to each other when they 
were nearest neighbors, lay within 30kms^^ of each other, and had logarithmic column densities 
differing by less than 1.0. 
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Fig. 8. — Top: RMS difference between the column densities of individually fitted absorption com- 
ponents in the two lines of sight versus the maximum measurement error permitted for the column 
densities in the sample used. Solid dots give the observed width of the A log N distribution, open 
diamonds the width expected on the basis of statistical errors in the column density determination 
alone. Bottom: the value of function E shows the trade-off between omitting measurement points 
with errors above a certain value, and the consequent reduction in the size of the sample, versus 
the maximum measurement error accepted in the sample. Positive values of E imply that it is 
beneficial to continue throwing away points with the largest errors, negative values imply a loss of 
information. The zero crossing corresponds to the optimum choice. 
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Fig. 9. — Distribution of differences between the HI optical depth At =| ta — tq \ for various 
column density ranges. The fraction of pixels in the column density ranges shown are 100% (solid 
line), 36.6% (dashed line), 36.4% (dash-dotted line), and 17.4% (dotted line). 



-34- 



wind goes off dense shell reaches LoS A 




LoS A 



LoSB 



pressure equilibrium LoS A gets disturbed 



spectrum A 



spectrum B 



if sound crossing time < shell crossing time: 



"v-^l/- 




Fig. 10. — left: A wind bubble expands near two lines of sight (A and B). The absorption spectra 
(bottom of the figure) look identical before the shell hits the first line of sight. Middle: the first 
line of sight (A) gets disturbed by the passage of the shell. The absorption pattern of A looks 
very different from that of B. Right: after the shell has slowed down sufficiently pressure waves 
re-establish equilibrium between the lines of sight. The absorption pattern although different from 
the initial state, looks again the same in both lines of sight. 
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Fig. 11. — Upper limit on the number density of galaxies surrounded by wind bubbles, as a 
function of the redshift zsf where the wind starts. The diagram is based on a 'disturbed' fraction 
IloS ^ 0.23 of the Lya forest. The different curves belong to different values of L^s/n^^. Winds 
at the low end of the L^s/n^^ range are less well constrained because they have a shorter lifetime 
for detection. The mean redshift of observation (3.26) is marked by an arrow. 
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Fig. 12. — Upper limit on the volume filling factor for wind blown bubbles. The input parameters 
are the same as for the previous figure. Strong winds L^^/n^rj ~ 100 — 1000 can fill a significant 
fraction of the volume of the universe (up to 40% if they start as early as redshift 10, up to ~ 18% 
if they start at z=4), whereas winds with L^s/n^^ < 1 are essentially unconstrained if they occur 
before redshift 7. 



